Pseudospin excitations in coaxial nanotubes. 
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In a 2DEG confined to two coaxial tubes the 'tube degree of freedom' can be described in terms 
of pseudospin- 1/2 dynamics. The presence of tunneling between the two tubes leads to a collective 
oscillation known as pseudospin resonance. We employ perturbation theory to examine the depen- 
dence of the frequency of this mode with respect to a coaxial magnetic field for the case of small 
intertube distances. Coulomb interaction leads to a shift of the resonance frequency and to a finite 
lifetime of the pseudospin excitations. The presence of the coaxial magnetic field gives rise to pro- 
nounced peaks in the shift of the resonance frequency. For large magnetic fields this shift vanishes 
due to the effects of Zeeman splitting. Finally, an expression for the linewidth of the resonance 
is derived. Numerical analysis of this expression suggests that the linewidth strongly depends on 
the coaxial magnetic field, which leads to several peaks of the linewidth as well as regions where 
damping is almost completely suppressed. 
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I. INTRODUCTION 

Due to the Coulomb interaction between the charge 
carriers, collective excitations called plasmons may result 
as the electromagnetic response in a solid-state structure. 
These excitations have been extensively investigated in 
bulk and low-dimensional systems from both theoretical 
and experimental points of view*i~— However, the inter- 
est on plasma excitations has been renewed due to the 
experimental possibility of tailoring interactions such as 
spin-orbit coupling, which may affect the properties of 
the plasmonsA"— The recent ability of producing novel 
two-dimensional systems such as graphene has also moti- 
vated new investigations of the plasmon dispersion.— ~— 

Interestingly, collective excitations can also emerge 
in systems involving two spatially separated two- 
dimensional electron gases which couple to each other 
through the Coulomb interaction^ - — A typical exam- 
ple is the excitation of different plasmon modes in bilayer 
systems, where two quasi-2D electron systems (each with 
only the lowest subband being occupied) are separated 
by a potential barrier. Even when the potential barrier 
is large and the tunneling is largely suppressed, the in- 
terlayer Coulomb interaction may couple the two quasi- 
2D electron systems if the interlayer separation (barrier 
width) is small enough. In such a case there are intralayer 
plasmon excitations in which the electrons in one of the 
layers may collectively oscillate in phase (optic plasmon 
mode) or out of phase (acoustic plasmon mode) with the 
electron oscillations in the neighboring layer . 19 ' 20 By de- 
creasing the barrier, the tunneling becomes relevant and 
splits the single subband in each quasi-2D electron sys- 
tem. As a result new excitations consisting of interlayer 
collective charge oscillations (intersubband or transverse 
plasmons) appear j^ r 18 i 21 Such interlayer collective oscil- 
lations have recently been re-interpreted as pseudospin 
excitations.— 

Within the pseudospin approach, the electrons in one 
of the layers are assigned one pseudospin state and the 
electrons in the other layer the opposite pseudospi m 22 ' 23 
The interlayer excitations can be regarded as pseudospin 



excitations mediated by the tunneling strength, which 
acts on the pseudospins as an effective magnetic field. 
Thus, analogous to the conventional ferromagnetic reso- 
nance in magnetized materials whose electron spins are 
manipulated by an external magnetic field, the tunneling 
(effective magnetic field) in the bilayer system leads to 
a pseudospin resonance describing the interlayer collec- 
tive mode (intersubband plasmon)^ -17 ' 21 Furthermore, 
since the pseudospin degree of freedom is an analog to 
the real spin, new pseudospintronic devices could be re- 
alized by means of controlled pseudospin manipulation, 
in close analogy with the control of real spin in spintron- 
ics applications . 2 ' 25 In particular, a pseudospintronic de- 
vice based on semiconductor bilayers has theoretically 
been suggested^ 2 , as the analog to the conventional spin- 
transfer oscillator i^r— 

Due to many-body effects the magnetization dynam- 
ics in magnetized systems are affected by the so-called 
Gilbert damping] 29 ' 30 Such effects also have their analog 
in bilayer systems, where the Coulomb electron-electron 
interaction produces a shifting of the pseudospin reso- 
nance frequency and leads to a finite lifetime of the ex- 
citations (that is, to damping). 22 Therefore, the investi- 
gation of the pseudospin excitations is of relevance for 
understanding the nature of correlations in bilayer-like 
systems. 

Another interesting issue is the investigation of the 
pseudospin resonance in systems with more exotic geome- 
tries. Nowadays techniques allow for the realization of a 
wide range of possible geometries by using semiconduc- 
tors or carbon based materials, for example. From this 
point of view, the phenomenon of pseudospin resonance 
offers the possibility of investigating many-body effects 
under different geometric configurations. In what follows 
we focus our discussion on the case of generic coaxial nan- 
otubes which can be experimentally realized from a vari- 
ety of materials including metals, metal-oxides, carbon, 
and semiconductors 4ir— Coaxial nanotubes are particu- 
larly interesting systems for pseudospintronics, since they 
exhibit both the bilayer-like behavior of pseudospin exci- 
tations and the interplay between many-body and coher- 
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FIG. 1: (color online). Schematic of two coaxial tubes in the 
presence of an external magnetic field B along their axis. The 
length of the tubes and the intertube distance are denoted by 
L and d, respectively. The radius of the inner (outer) tube is 
given by R - d/2 {R + d/2). 



cnt effects which become apparent when a coaxial mag- 
netic flux pierces the system. In such a case, as will be 
shown below, pseudospin resonance can be induced not 
only by an external electric field but also by fluctuations 
of the coaxial magnetic field. 

Although there have been some investigations on 
plasma excitations in coaxial nanotube a 39 ' 40 these studies 
were limited to the non-tunneling regime treated within 
the random-phase approximation (RPA). Here, we con- 
sider the possibility of tunneling (and therefore of pseu- 
dospin excitations) between the inner and outer tubes. 
Furthermore, we use a perturbative scheme recently pro- 
posed in Ref. [22| which appears to be superior to the 
RPA. 

The paper is organized as follows: In Sec. HU the pseu- 
dospin degree of freedom and the model Hamiltonian of 
the system are introduced. To construct the Hamilto- 
nian, we start with a single-particle Hamiltonian tak- 
ing into account tunneling effects and analyze the en- 
ergy spectrum of the resulting, before we include the 
electron-electron interaction. In Sec. IIII1 the pseudospin 
resonance and its corresponding response function are in- 
troduced. The perturbation scheme is set up in Sec. IIVI 
The scheme is then applied in Secs.lVland lVII to calculate 
the resonance frequency and the corresponding lincwidth, 
respectively. The paper is concluded by a small summary. 



II. THEORETICAL MODEL 

The system (see Fig. [IJ consists of two coaxial tubes, 
which have the radii R±d/2 and the length L 3> R 3> d. 
A 2DEG is confined to the surface of each cylinder. 

To approximate the band structure of actual solids, 
we assign an effective mass m to the electrons. Similarly, 
the Coulomb interaction contains the background dielec- 
tric constant e, which is defined by e = e r eoi where e r is 
the relative dielectric constant of the solid. The corre- 
sponding Bohr radius in this material is then defined as 
as = 4:ireh 2 /me 2 . The interior of the cylindrical system 



is threaded by a static, coaxial magnetic field B, which 
acts as a control parameter. To describe paramagnetic 
effects induced by the magnetic field, that is, Zeeman 
splitting, the g factor of the material has to be intro- 
duced. In this model electrons can either be located on 
the outer or inner tube. 



A. Single-particle approximation 

In the absence of tunneling (tunneling effects will be 
included later on) the wave functions are localized at the 
tubes and can be approximated as 



*out/in (r) = $out/in (f, z) /out/in (p), 



(1) 



where we have used cylindrical coordinates (see Fig[l}. 
The radial localization is characterized by the function 



out/in 



(p)\ 



s(p-rtI) 



(2) 



where S(x) represents the Dirac-delta function and the 
radial quantum number denotes whether the electron is 
located on the outer (out) or inner (in) tube. The lon- 
gitudinal and azimuthal motions as well as the physi- 
cal spin of the electrons are described by the spinors 
3> ou t/i n {<p, z). Apart from the physical spin, the radial 
motion can be reinterpreted as an effective two-level sys- 
tem, which we describe in terms of pseudospin- 1/2 dy- 
namics. We then replace the functions / ou t/in by pseu- 
dospinors whose labels ft and ft correspond to the wave 
function localization in the outer and inner tubes, respec- 
tively. 

In the absence of tunneling the pseudospin system is 
described by the single-particle Hamiltonian 



Ha 



(3) 



where 



(ft* 



* - a yu, r ) is the operator of annihilation 

7 7 i Th . rC . (7 1 A 

(creation) of a particle with angular momentum n, mo- 
mentum along the z-axis k, spin a =t, 4- aud pseudospin 
A =ft, ft. The single-particle eigenenergies are given by 
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Here, and in what follows, we use capital and lower- 
case Greek characters for denoting pseudospin (ft, JJ.) and 
physical spin (t,|), respectively. We have also intro- 
duced the average magnetic flux $ = ttBR 2 , the fluxon 
$0 = h/e, and the free electron mass m e . 

We now consider the possibility of uniform tunneling 
between both tubes. The tunneling amplitude A is as- 
sumed to be independent of the external magnetic field. 
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Such an approximation is reasonable in systems in which 
the confinement is stronger than the cyclotron effects. 
The tunneling Hamiltonian, 



Hi 



(5) 



with a x as the corresponding pseudospin Pauli matrix, al- 
lows for the coupling of states localized in different tubes 
but with the same values of the quantum numbers a and 
n, which is consistent with the conservation of spin and 
angular momentum during the tunneling. 

The form of Ht [see Eq. ([5])] makes it clear that the 
tunneling amplitude can be interpreted as the pseudospin 
analog to a magnetic field in the ^-direction. This term 
arises due to the overlap of the actual radial wave func- 
tions (which in reality are not as perfectly localized as 
our (5-like model functions). 

In analogy to the spin operator we can introduce the 
pseudospin vector operator S whose ^-component, 



characterizes the tunneling between the tubes. Indeed, 
the tunneling Hamiltonian in Eq. (J5J) can be rewritten in 
second quantization as 

H t = -AS x ,. (7) 

The y-component of the pseudospin operator is given by 



(8) 



and measures the tunneling current flowing between the 
two tubes, while the z-component, 



2 / j y n,k,o 



(9) 



measures the charge accumulation in the tubes. In 
Eq. ([SJ N-ft and denote the number-of-particles op- 
erators in the outer (ff) and inner (JJ-) tubes, respectively. 

The single-particle Hamiltonian including the tunnel- 
ing effects is given by Hq + H t . Since S z and S x do not 
commute, the quantum number A =ff-, J|, which classi- 
fies the eigenvectors of S z , has to be replaced by another 
quantum number A = ±1 describing the pseudospin de- 
gree of freedom in the presence of tunneling, whereas 
the momentum along the z-axis as well as the angular 
momentum and spin projections on the z-axis remain 
good quantum numbers. The new pseudospinors with 
pseudospin quantum number A are linear combinations 
of the pseudospinors with A =f|~, JJ, and describe bonding 
(A = —I) and antibonding (A = +1) states. In the basis 
of A-pscudospinors, the Hamiltonian Hq + H t becomes 
diagonal with the energy spectrum given by 



£ra,fc,<r,A 



e n,k,a,it + e n,k,<r,ty 



A 



-A r , 



(10) 




FIG. 2: (color online) Single-particle energy spectrum in the 
presence of tunneling for different values of the magnetic 
flux [(a): $ = 0, (b): $ = 0.05$o, (c): $ = 0.10$o, (d): 
$ = 0.15*o, (e): $ = 0.20$o, (f): $ = 0.25$o] with parame- 
ters A = 0.1eF, re f , d = 0.1-R, g = 2, r s = 1, R = 0.5asr s , and 
m = m e . The quantum numbers n, a =t,4-i an d A = ±1 re- 
fer to the angular momentum, physical spin, and pseudospin, 
respectively. In all the cases t n ,h,a,\ > £n,k,a,-i, as explicitly 
indicated in (a). 



where 



A„ = \J (e^fc^jj, - £ n ,k,<j^) 2 + A 2 . (11) 



Here, e n ,k,<T,it an d en,fc,<7,4 are the eigenenergies of Hq 
which are given by Eq. (j4]). 

From Eq. (1101) one can see that the energy spectrum 
consists of one-dimensional subbands labelled by the 
quantum numbers n, a, and A. The evolution of the 
lowest subbands with increasing magnetic flux is shown 
in Fig. [2] The quantity of reference in this plot is eF,rcf , 
the Fermi energy of a flat 2DEG with the same den- 
sity n e as the cylindrical 2DEG considered here, and 
&F,ref = ^2meF,rcf/fr 2 - We assumed a model system 
with parameters A = O.leF.rcf, d = 0.1R, g = 2, r s = 1, 
R = 0.5asr s , and m — m e . Here, r s = \/ yJnn e a 2 B is the 
Wigner-Seitz density parameter. Note, that in Fig.[2]thc 
energy spectrum is only shown in the zeroth order in d, 
because in the final results of our calculations the single- 
particle eigenenergies enter only in the zeroth order. 

At zero magnetic field [see Fig. Ufa)] each subband - 
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except the subbands with n = 0, which are doubly de- 
generate - has a fourfold degeneracy due to spin degen- 
eracy and rotational invariance. For the chosen param- 
eters, two degenerate pairs of subbands, namely a pair 
denoted by n = 0, a =4-, A = ±1 and a pair denoted by 
n = 0, a =t, A = ±1, are occupied as shown in Fig. [2ja). 

A finite magnetic field lifts the spin as well as the angu- 
lar momentum degeneracy. With increasing strength of 
the magnetic field the energies of the two subbands with 
a =t start to increase towards the Fermi level, while 
the energies of the two subbands denoted by a =4- de- 
crease. This results in the occupation of four, five, and 
six non-degenerate bands, as shown in Figs. [2fb), (c), 
and (d) , respectively. Further increasing of the magnetic 
field strength leads to the inversion of the order of the 
(n = 0, a =t) and (n = —1, a =4-) subbands [compare 
Figs. Hfd) and (e)] and, eventually, to the depopulation 
of the (n = 0, a =t) subbands [see Figs.[5ff)]. 



B. Effects of Coulomb interaction 

We now include the effects of the electron-electron in- 
teraction. Electrons in the same tube interact via the 
intratubc potentials 



^o U t/i„ (hq) = -Ii \q\R± kb Ki \q\R±\q 



'2 J ' 



(12) 

depending on the tube both electrons are located in. On 
the other hand, the intertube Coulomb interaction be- 
tween electrons from two different tubes is given by 

V d (I, q) = (\q\R - \q\fj K t (\q\R + \q\fj . (13) 

In these expressions, h{x) and Ki(x) are the modified 
Bessel functions, while I and q denote the change of angu- 
lar and linear momentum (along the z-axis), respectively. 
It is convenient to define linear combinations 

V t t (I, q) = ~ [V out (I, q) ± V d (I, q)] , (14) 



and 



V±(l,q) = -{V in (l,q)±V d (l,q)}. (15) 



Then, the electron-electron interaction can be written as 



H' = 

E { Kut (I, 9) + Vr (I, q)] S z (I, q) S z (-/, -q) 

+ I fc t (I, q) + V+ (I, q)} h (I, q) h (-/, -q) 
+ 3 Kut (I q) ~ Vr (I, q)] [S, (I, q) h (-1, -q) 
+n(l,q)S z (-l,-q)]} 

-dzEfe (i,q)-vr (l,q)]S z , 

l,q 



where 

n(l,q) = a n,fe,a,A a n+9^+'^.A (17) 

is the local density operator. The complete Hamiltonian 
describing our system is 



H = H + H t + H', 



(18) 



which comprises single-electron, tunneling, as well as 
Coulomb coupling terms. 



III. PSEUDOSPIN RESONANCE 

Our goal is to investigate the pseudospin resonance, 
which is an analog to the ferromagnetic resonance ) 29 i 30 
and how the resonance is affected by the electron-electron 
interaction. In what follows, zero temperature is con- 
sidered. In the uniform case, that is, for zero momen- 
tum and angular momentum transfer, the pseudospin re- 
sponse function is, in the linear response regime, 



1 



2ttL 



"2 j S 2 



(19) 



where the brackets denote the Kubo product, 



((A,B)) u = -i / dr 



[A(t),B(0)]\0) 



(20) 

Since S z measures the difference between the number of 
electrons in the outer and inner tubes, the pseudospin 
resonance describes collective oscillations of the particle 
densities between the tubes. This collective mode can be 
induced either by an external electric potential V ex t (t) 
applied between the tubes (as is also the case in flat- 
bilayer systems) or by fluctuations Bfi (t) in the coaxial 
magnetic field (or by applying an oscillating, coaxial aux- 
iliary magnetic field) . The latter case is dealt with in the 
model by replacing the constant magnetic field ampli- 
tude B with B + Bfi(t) and treating the arising linear 
term containing Bji(t) as an external perturbation and 
neglecting the higher order terms. Using linear response 
theory-'* the fluctuations of the pseudospin expectation 
value of the system due to those perturbations can be 
calculated from the following expressions: 



5 z (u)) = ~ X (u)Vest(M>) (21) 



and 



6&( U )) = ~±x(t»)B,d«,), (22) 

where T4 x t (w) and Bfi{u) are the Fourier transforms of 
V cxt (t) and Bfi(t). By comparing Eqs. (|2"Tj) and (|2"2")l one 
can see that 



(16) 



V B (u>) 



h d $ 

2mfl¥^ 



(23) 
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acts as an effective field which has the same effect as an 
external electric potential. 

The pseudospin resonance is given by the condition 



Re [x^ 1 Kcs)] = 0, 



(24) 



from which the resonance frequency w res will be ex- 
tracted. 

As suggested by Eq. (|2"21 , in coaxial tubes an auxiliary 
magnetic flux piercing the system leads to fluctuations 
of the pseudospin expectation value and can therefore 
be used as an alternative control parameter for inducing 
collective oscillations. 



IV. PERTURBATION THEORY 

We will work in the limit d <C R and introduce the di- 
mensionless intertube distance ^ =d/flasan expansion 
parameter. To calculate the resonance frequency in pow- 
ers of £d, a slightly modified version of the perturbation 
theory developed in Ref. [H is applied. The perturbation 
scheme from Ref. [22| has been modified to account for 
the cylindrical geometry considered here. First, x ( w ) is 
systematically decomposed into ground state expectation 
values and correlation functions, and an exact equation 
for x analogous to Eq. (3) in Ref. [22l. is derived. This 
scheme is rather elaborate and we refer to the Appendix 
lAl [see Eq. (|A3|) ] for more details. The perturbation the- 
ory is based on the expansions 



KTut (I v) = e -^i'i (\i\R) k x (\ g \R) d + o(£) 

e 



and 



Vr (l,q) 



-^I l (\q\R)K' l Qq\R)d+0(&), 



(25) 



(26) 



which vanish in the zeroth order and can be considered 
small perturbations for ^ = d/R <C 1. To calculate 
the correlation functions and expectation values, it is 
convenient to switch from the Heisenberg to the inter- 
action picture and base the perturbation scheme on the 
perturbation H per , which consists of those parts of H' 
that contain at least one factor V^~ t (l,q) or V- m (l,q). 
The unperturbed Hamiltonian -ff U npor is then given by 
H, with V^~ t (l,q) and VJ~ (l,q) set to zero. In the fol- 
lowing, we will restrict ourselves to the high-density limit, 
where we can neglect the density-density coupling in H' 
(which is not part of H pcr ) and use the Hamiltonian 

-Hunper = Hq + H t and its ground state as our reference 
system. 

The strategy is to expand \ (uj) in powers of £d (see 
Appendix lAl for details). The expansion series of xfa) 
is then inverted and expanded again in powers of £d up 
to the same order as x ( w ) was - After that, the zeros of 
this new power series are computed as a function of the 
intertube distance and expanded once again , which pro- 
vides the expansion series of the resonance frequency in 



powers of As an important example, the pseudospin 
resonance frequency will be calculated in the following 
section. 



V. PSEUDOSPIN RESONANCE FREQUENCY 

An approximate expression for the resonance frequency 
can be obtained by expanding the pseudospin response 
function up to the first order in In such a case the 
evaluation of x (w) [see Eq. (|A3|) ] is greatly simplified (see 
Appendix [A] for details) and one obtains the following 
approximate relation, 



2A 2 



xEfc('.9)+^(i.«)] w / (W (i.«) (27) 

l,q 

t2s 



where 



and 



M x = 



2ttL 



<0| S x |0) 



(28) 



(29) 



Expanding j\A x in powers of £d one obtains in the zeroth 
order 



M 0) 



irT, Z-^ 



,(0) 



,k.<7, — 1 



1 (0) 

L n,k,a, + 1 



2ttL 



with 



,(0) _ ft (0) (0) 



(30) 



(31) 



where uip^ = e^' / h [cp 1 is the Fermi energy of the sys- 



J0) 



tern in the zeroth order in ^] , and ui^ — ■ 



,fc,CT,A 



n,fc,(7,A 



/h. 

The correction Ai x ^ adding the first order terms of M. x 
is irrelevant for the expansion of x _1 (uj) up to the first 
order in £<j- I n the equations above, and in what fol- 
lows, the superscripts denote the respective order in 
Finally, we have introduced 



L n+l,k+q,a,A 



/<°> (l,q) =6i, 6 g , [2*LM<® 

+ 1 V n (0) (l- 

n,k,a, A 



l {0) 

l n+l,k+q,a-,'S: 



n.k.c A,S 



in Eq. (gT}. 



(32) 
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Next, the approximate pseudospin response function 
given in Eq. (1271) is inverted and then expanded in pow- 
ers of £d- From the zero of this series the resonance fre- 
quency uj rcs can be calculated. After several algebraic 
manipulations, the resonance frequency reads as 



A 2 



K 2 h 2 e R 



4Ae 2 dM x 0) (t)MX } (I) 



(0) 



M 



(«) 



0{e d ). (33) 



For the derivation of this expression, the Wronskian for- 
mula 



Ii(x)Kl{x) -Il^K^x) 



1 



(34) 



for modified Bessel functions^ has been used to rewrite 
the potential differences in the first order in 

fc(U) + ^(U)] (1) = y| (35) 

In Eq. (|33|) we have introduced the spin-resolved pseudo- 
magnetizations M x °^ (t) an d M x °^ (I), which are defined 

as 



x v ; 2ttL ^ 



l n,k,a,-l 



l (0) 

L n,k,a, + 1 



(36) 



i ; , k 



Recall that the pseudospin quantum number A = ±1 
refers to the bonding and antibonding eigenstates of the 
Hamiltonian Hq + H t with eigenenergies t n ,k.a,h, given 
by Eq. (fTU)) . In Eqs. (l3"Tj) and (j3"6"j) the energy spectrum 
enters only in its zeroth order in f^. 

It follows from Eqs. (|31l) and (f31)f that the spin- 
resolved pseudomagnetization measures the difference 
between the number of bonding (A = — 1) and antibond- 
ing (A = +1) states for a fixed spin quantum number a. 
Therefore, the total pseudomagnetization is given by the 
sum 



M 



(o) _ 



E 

r =tA 



(37) 



From Eq. (|33l) one can see that the Coulomb inter- 
action produces a shift of the resonance frequency with 
respect to A/h, the resonance frequency one would have 
obtained within the single-particle picture. This opens 
the possibility of investigating many-body effects by mea- 
suring the frequency shift of the pseudospin resonance. 

At finite magnetic fields, the shift in the resonance fre- 
quency, given by 



A 2 

h 2 



(38) 



depends on the different spin populations of the subbands 
denoted by the pseudospin quantum numbers A = ±1. 
This is due to Zeeman splitting which lifts the spin degen- 
eracy of the different subbands. The division into spin- 
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FIG. 3: Magnetic field dependence of the shift in the reso- 
nance frequency for a model system with parameters A = 
0.1eF, rc f, d = 0.17?, g = 2, r s = 1, R = 0.5asr s , and m = m e . 



important if the Zeeman term is not small compared to 
A. 

The magnetic field dependence of the shifted resonance 
frequency £1 is shown in Fig. |3] for a model system with 
parameters A = O.lep.rcf, d — 0.1R, g = 2, r s = 1, 
R = 0.5asr s , and m = m e , that is, the same parame- 
ters as in Fig. [2] One can see that for small values of 
the magnetic field the shift of the resonance frequency is 
almost constant. This trend changes abruptly when the 
magnetic flux reaches the value 0.086 $o- At this point 
of non-analyticity the resonance frequency shift sharply 
starts to increase with the flux. For larger magnetic fields 
a pronounced peak of f2 2 develops. However, increasing 
the magnetic field even further results in a sharp drop of 
f2 2 , and the shift of the resonance frequency vanishes for 
higher fields. 

The behavior of the resonance frequency shift is de- 
termined, essentially, by the magnetic field dependence 
of the spin-resolved pseudomagnetizations [see Eq. (|3"3"|) ]. 
which are shown in Fig. |4j As already mentioned, 
the spin-resolved pseudomagnetizations (t) an d 

Ai^ (4-) measure the difference between the number of 
occupied bonding (A = —1) and antibonding (A = +1) 
states for up-spin and down-spin particles, respectively. 
Therefore, the magnetic field dependence of the spin- 
resolved pseudomagnetization can be qualitatively ex- 
plained by analyzing the changes of the energy spectrum 
with the magnetic flux [see Fig. [2j. 

In the limit of zero magnetic field both the bonding 
and antibonding states are spin degenerate [see Fig.[2ja)]. 
Therefore, the single-particle energies no longer depend 
on the spin, implying 



M°>(t) = M 0) a) = ^M 0) , ( 39 ) 



resolved pseudomagnetizations (t) an d -Mx'* (4-) is which can be seen in Fig.0J In such a case the resonance 
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FIG. 4: (color on line). Magnetic field dependence of the spin- 
resolved and total pseudomagnetizations for the same generic 
system as in Fig. [3j 

frequency shift given by Eq. (|3"3"j) reduces to 

& = ^±MW+0(e d ). (40) 

Note that since the eigenstates are non-degenerate with 
respect to the quantum number A [Fig. Ufa)], the bond- 
ing and antibonding states are unequally occupied. This 
results in a finite value for the spin-resolved and total 
pseudomagnetizations at zero magnetic field. 

The relation in Eq. (|^U|l is similar to that of the flat 
bilayer system examined in Ref. |22L In fact, if we com- 
pare the different constituents to the frequency shift in 
Eq. (|40|) and in Ref. and their respective definitions, 
we can see that, in the first order in £d, the structure of 
both expressions is the same. 

At small magnetic fields the energy shift (and, there- 
fore, the difference in occupation) of the bonding and 
antibonding states for both up- and down-spin particles 
changes smoothly with the field strength, leading only to 
small effects on the spin-resolved pseudomagnetizations. 
Nevertheless, one can still observe an increase (decrease) 

in M { x } (t) [A4 0) (I)] as the non-degenerate subbands 
with n = 0,a- =f,A = ±1 (n = 0,a =|,A = ±1) shift 
up (down) towards (from) the Fermi level and the dif- 
ference between their populations increases (decreases). 

The opposite behavior of Mx (t) and (I) com- 

pensate each other, resulting in almost magnetic field- 
independent total pseudomagnetization and resonance 
frequency shift, in the region < $/$o < 0.086 [see 
Figs. [3] and 2]. The situation changes drastically at 
$ ~ O.OSB^o, at which point a fifth subband, namely 
the subband with n = — 1,(7 =4, A = —I, becomes 
occupied [see Fig. [2jc)]. This results in a sharp in- 
crease of (I) because now there is an additional 
subband, which is rapidly populated and contributes to 
the A = — 1 states in M x °' (4), while there is only one 
band with A = +1 that contributes to (l). At these 



FIG. 5: Magnetic field dependence of the shift in the reso- 
nance frequency for a narrow model carbon nanotube system 
(e r = 2.4, m = 0.25m e , g = 2, R in = 1.1 nm, R out = 1.44 nm, 
ni = 9.1 nm -1 , and A = 0.1eF, rc f) 



magnetic fields, the spacing in the energy spectrum be- 
tween states which differ only in their pseudospin quan- 
tum number A is already much smaller than the spac- 
ing between states which differ only in their real spin 
quantum number a. At $ « 0.117$0j the subband with 
n = —1,<7 =|, A = +1 becomes occupied [see Fig. |2fd)], 
which consequently leads to a decay in the absolute value 
of A4^ (I)- This behavior continues until the subbands 
with n = —1,(7 =i,A = — I and n = —l,a =4, A = +1, 
which are energetically very close to each other, lie far be- 
low the Fermi level [see Fig.[2jf)], and the decay becomes 

less pronounced. If we look at the behavior of Ai^ (t) in 
Fig. 21 we can observe a steady increase, which becomes 
more steep as the minimum of the antibonding subband 
with n — 0, a =t, A = +1 moves closer to the Fermi level 
[see Figs. HJc) and (d)]. At $ w 0.21<J> this subband 
is no longer occupied, while the bonding subband with 
n = 0,a =f,A = -I is still populated [see Fig. HJe)]. 
This results in the peak of (t)> which can be seen 
clearly in Fig. [4] With even higher magnetic fields the 
remaining bonding subband with n = 0,(7 =t 5 A = — I 
moves up and becomes less and less populated, leading to 
the decay of (t) until the band minimum crosses the 
Fermi level (at $ w 0.234$ ). For $ > 0.234$ the spin- 
up subbands lie already above the Fermi level and only 
spin-down subbands are occupied [see Fig. HJf)]. Conse- 
quently, (f) vanishes in this region [see Fig. [3]. 

Finally, the magnetic field dependence of the resonance 
frequency shift [see Fig. [3] is determined by the interplay 
of the above discussed behaviors of both the spin-resolved 
and total pseudomagnetizations. 

In Figs. [5] and [5]we show the magnetic field dependence 
of the resonance frequency shift for realistic model sys- 
tems and for magnetic fluxes which correspond to mag- 
netic fields of up to 40 T. In Fig. [5] the system parame- 
ters e r — 2A, m = 0.25m e and 5 = 2 have been chosen 
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FIG. 6: Magnetic field dependence of the shift in the reso- 
nance frequency for an InGaAs/GaAs nanotube system (ob = 
9.8 nm, g = -0.44, R = 10 nm, d = 1 nm, n = 1 x 10 11 cm -2 , 
and A = 0.1e F rct ) 



to simulate a narrow carbon nanotube system with the 
radii i?,j n = 1.1 nm and i? ou t = 1.44 nm and the electron 
density per unit length m = 9.1 nm -1 ! 33 ' 39 ' 42 " — 

In Fig. |5] the parameters correspond to a model In- 
GaAs/GaAs nanotube system (as = 9.8 nm and g = 
—0.44) with an average radius of R = 10 nm, a distance 
d = 1 nm between the two tubes, and an electron density 
n = 10 11 cm- 2 ^M 

As in Figs. 131 and |H the magnetic field dependence of 
the pseudomagnetizations and resonance frequency shift 
displayed in Figs. [5] and [6] can be explained by analyz- 
ing how the single-particle energy spectrum evolves with 
the external magnetic field. The only difference is that 
in Figs. [5] and [5] more subbands are occupied than in 
Fig. [3l which results in a more complicated structure of 
il 2 . Furthermore, the effect of Zeeman splitting is not 
as strong as in Fig. [3J where the parameters were chosen 
for simplicity to illustrate Eq. [33] Therefore, one would 



have to go to unrealistically high magnetic fields to ob- 
serve the vanishing of the resonance frequency shift. For 
magnetic fields below 40 T one has not yet entered this 
regime for the systems shown in Figs. [5] and [5] There 
are several peaks for magnetic fields below 40 T in the 
InGaAs/GaAs nanotube system (see Fig. [5]). In the nar- 
row carbon nanotube system, on the other hand, the res- 
onance frequency shift does not vary much for the fields 
considered (see Fig. [5]) . Only at unrealistically high mag- 
netic fields pronounced peaks would appear in this sys- 
tem and one would need even higher magnetic fields to 
observe the vanishing of the shift of the resonance fre- 
quency. 

VI. LINEWIDTH OF THE PSEUDOSPIN 
RESONANCE 



Additional corrections to Eq. (|3"3")l are obtained by ex- 
panding the pseudospin response function up to the sec- 
ond order in at which point a finite imaginary part of 
X^ 1 (w) emerges. A non- vanishing imaginary part implies 
damping and a finite lifetime of the mode. This means 
that the pole of x ( w ) is -at least up to the second order in 
£<i- replaced by a roughly Lorentzian-shaped peak at the 
resonance frequency. The linewidth T of the resonance 
can be obtained from the imaginary part of the zero of 
X~ x {ui) (in contrast to the zero of Re [x^ 1 ( w )] > which 
only yields the resonance frequency). This situation is 
the pseudospin analog to the finite linewidth of the fer- 
romagnetic resonance induced by the Gilbert damping. 

Using the shorthand notations u = (n, k), v = (l,q), 
and v' = (l',q'), the linewidth up to second order in ^ 
is given by 



3 2 d 



lim 



A{co) B{uj) 



Here, A (w) and B (oj) denote the sums 



(41) 



^h = i^EEE^-^ (i - nffi™*) (i - »£U»S^,-.,-a (42) 

v / u,cr v,v' A,S 

and 



B{Uj) ^2l£y £ E V 2 (v)^-^ ( v,v')) (l-n£U,,A) (l-41 W Kt, 

^ ' U,V,V 7 (7,(7' , A 

£ Ev(v)v ( v')^-^(v,v')) (i-»£Ua) (i-4 2w 

> ' 11. V. v' tT.A 



(0) 

A"-u+v-v',<j',A 



X] n u,cr,A n u+v-v',o-,A i 

(43) 



where ~- " h ... h 



^A(v,vO = ^«r+A^ + ^ A , (44) 
ZmK m n 



n ^ = ^U> + -M', (45) 
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FIG. 7: Magnetic field dependence of the linewidth of the 
pseudospin resonance for an InGaAs/GaAs nanotube system 
(an = 9.8 nm, g = —0.44, R — 10 nm, d — 1 nm, L = 1 fim, 
n = 1 x 10 11 cm" 2 , and A = 0.1e%) 



and 

V (v) - | 9 | [// K t (\q\R) + I, (\q\R) K[ (\q\R)] . 

(46) 

The zeroth order band occupation numbers n ® « have 
already been introduced in Eq. pip . The physical origin 
of the damping can be understood by looking at the ex- 
pressions for A (ui) and B (w). The mode loses energy by 
exciting two particle-hole pairs out of the single-particle 
spectrum. The process described by A (uj) consists of an 
excitation where the pseudospin is conserved and another 
one where the pseudospin is nipped. Furthermore, the 
total momentum, the total angular momentum and the 
total spin are conserved during this double-excitation. 
The second process is given by B(tv), which also de- 
scribes a double-particle-hole excitation. Like before, the 
process is momentum-, angular momentum-, and spin- 
conserving. But in this case the process involves cither 
only bonding or only antibonding states, and there is no 
pseudospin-nip. The process given by A (lu) can also be 
found in the flat bilayer system^ while B (lu) is a man- 
ifestation of the cylindrical system and the existence of 
two different intratube-potentials. 

The magnetic field dependence of the resonance de- 
cay rate computed from Eq. (|41[) is shown in Fig. [7] for 
the case of a semiconductor system with the same pa- 
rameters as in Fig. [5] Additionally, we chose a finite 
length L = 1 fxm for the nanotubes to speed up the 
time-consuming numerical evaluation of Eq. (|41[) by dis- 
cretizing the linear momenta in Eq. (|41[) . The decay 
rate exhibits a strong dependence on the magnetic field, 
with several pronounced peaks at which the resonance 
linewidth is enhanced and regions where the damping 
is almost suppressed (that is, where T is close to zero) 
and the pseudospin resonance becomes very sharp. This 
interesting behavior makes the external magnetic field, 
which is an experimentally tunable parameter, attrac- 



tive for the controlled switching of the damping of the 
pseudospin oscillations. 



VII. CONCLUSION 

We have considered a cylindrical bilayer system con- 
sisting of two coaxial tubes. To account for the presence 
of two tubes, we have introduced a new quantum number 
which describes the two-level system. This two-level sys- 
tem was then interpreted in terms of pseudospin- 1/2 dy- 
namics. We have incorporated tunneling between the two 
tubes and calculated the energy spectrum of the single- 
particle Hamiltonian. The inclusion of tunneling be- 
tween the two tubes made it possible to find a pseudospin 
analog to the ferromagnetic resonance. Taking into ac- 
count the electron-electron interaction, the pseudospin 
resonance frequency has been calculated up to the first 
order in the intertube distance by using a perturbativc 
scheme, which has been developed for a flat bilayer sys- 
tem and appears to be superior to the RPAi22 Due to the 
Coulomb interaction there is a shift in the resonance fre- 
quency, which is also dependent on the coaxial magnetic 
field. This dependence results in pronounced peaks of 
the resonance frequency shift at certain magnetic fields. 
The shift, however, disappears at higher magnetic fields. 
The damping effects induced by the Coulomb electron- 
electron interaction on the pseudospin resonance have 
been investigated by computing the linewidth of the res- 
onance. The linewidth exhibits a strong dependence on 
the magnetic field, with a multi-peak structure. Apart 
from the peaks, where the decay rate is enhanced, in some 
ranges of the magnetic filed strength the damping of the 
pseudospin oscillations is almost suppressed. 
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Appendix A 

Here, we present the expression for \ ( w ) 011 which the 
perturbation scheme is based. The application of the 
Kubo product formulas 

((A,b)) =i<0|[A(0),B(0)]|0> + l((ftA,B)) , 

(Al) 

and 

(( i ^». = ^ l[ i ( )'^ )]l°)-^(( i '^».' 

(A2) 
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which can be verified by partial integration, leads to the 
following expression for the pseudospin response function 

xM = 



A iA 

M x - - r „ M d (w) 



m 2 2TrLh 2 ni 

A 1 





(2irL) 2 


iAui 


1 


h 2 nt ( 27r L) 2 


A 2 


1 




(2ttL) 2 


A 2 


1 




(2ttL) 3 


A 2 


1 




(2ttL) 3 


A 2 


1 




(2ttL) 3 


A 2 


1 


( 27r L) 3 


iA 2 


1 




(2ttL) 2 



2A 



^V 1 (v)(^.g(v) + 1 -/(v) 

V ^ 

5>2(v) (/i„ (v) — 2wLAiy) 

V 

^V 2 (v) (ft, (v)-2iLM 2 ) 

V 

$3EVi(v)V 1 (v')£o (vy.w) 

V v' 

$3EVi(v)V 2 (v')A (v,v» 

v v' 

V v' 

^5> 2 (v)V 2 (v')A (vy,«) 

V v' 

(Vi(v)Ca (v, w) + V 2 (v)C 3 (v, w)) , 

(A3) 

where we have introduced the quantities 

Vi(v)= [Kut (v) + (v)] , (A4) 



V 2 (v)= [V- t (v)-Vr ( v )], 

I Si |0) 



A4, 



27!\L 



t=o 



(A5) 
(A6) 



/ (v) = (0| S x (v) S x (-v) - 5 Z (v) 5 Z (-v) |0) 
g (v) = (0| S y (v) ^ (-v) + S z (v) 5„ (-v) |0> 



hi (v) = - (0| Si (v) n (-v) + h (v) $ (-v) |0) 



t=o 



C (v,v',w) = US xz (v),S xz (v' 
A (v,v» = ((i 2 (v) I Q(v') 
£ (v,v',o;) = /(q(v),Q(v') 



(A7) 

t=0 

(A8) 

(A9) 
(A10) 

(AH) 

(A12) 



A (v,v',a;) = ((Q(v) ) 4 z (v / ))) w , 
A (w) = 



(A13) 
(A14) 
(A15) 
(A16) 

S xz (v) = S x (v) S z (-v) + 5 2 (v) S x (-v) , (A17) 

Q(v) 



C 2 (v,w) = ^4«(v), [$,,#0])) 
C 3 (v,u>) = ((Q(v),\s y ,H 



l r 

2 



2S X 



S x (v) n (-v) + n (v) (-v) 

(A18) 

and the shorthand notations v = (I, q) and v' = (l',q'). 
This equation is analogous to Eq. (3) in Ref. but due 
to the intrinsic dependence of the single-particle eigen- 
states of Ho on the interlaycr difference and different in- 
tcrlayer Coulomb potentials in each tube, there are addi- 
tional terms. If this intrinsic interlayer dependence were 
absent, both intralayer Coulomb interactions would be 
identical, which means every term that contains the fac- 
tor [V^ (I, q) — (I, q)] would vanish. Furthermore, 

Ho would no longer depend on d and would commute 
with S y. I n this case, Eq. (|A3|) would reduce to Eq. (3) 
in Ref.lH 

To obtain the resonance frequency in the first order, 
we need to expand the response function in Eq. (|A3|) up 
to first order in In such a case, we do not have to 
consider the terms containing the products Vj(v)Vj(v') 

j = 1,2) because they are already of second order or 
higher in £<j. The situation is further simplified by look- 
ing at the commutator of Ho and S y , which also vanishes 
in the zeroth order in ^ and thus can only contribute to 
first or higher order of the pseudospin response function. 
This implies, however, that products of this commutator 
with Vi(v) (i — 1,2) do not add to the pseudospin re- 
sponse function up to the first order. As a result, many 
terms in Eq. (|A3|) can be omitted in the first order. The 
calculation of the remaining terms shows that only f(l, q) 
and the pseudomagnetization A4 X do not vanish in the 



Sy, HQ 



■ S z 



also 



zeroth order in Moreover, 

vanishes in zeroth and first orders in which means 
that this term cannot contribute to the pseudospin re- 
sponse function up to the first order. Consequently, we 
are left with the formula given in Eq. (|33[) . which is a 
simplified expression of the response function accurate 
up to the first order in £<j. 

For the linewidth of the pseudospin resonance, one has 
to calculate x ( w ) m a t least second order in £<ji which 
means no term in Eq. (|A3[) can be omitted. 
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